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Abstract 

A finite expansion of the exponential map for a, N x N matrix is presented. The method 
uses the Cayley-Hamihon theorem for writing the higher matrix powers in terms of the first 
N-1 ones. The resulting sums over the corresponding coefficients are rational functions of the 
eigenvalues of the matrix. 
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1 Introduction 



In the Lie theory of groups and their corresponding algebras the exponential map is a crucial 
tool because it gives the connection between a Lie algebra element € g and the corresponding 
Lie group element g G G 



exp 




(For details see Q and |T|, |2[ and ref. therein). In some low dimensional cases, like SU(2) and S0(3) 
the explicit expansion of the exponential map is known. Some years ago the exponential map 
for the Lorentz group was given by W. Rodrigues and J. Zeni |8|. For the higher dimensional 
groups SU {2,2) and 0(2,4) a method for the expansion was developed by A.O. Barut, J. Zeni, 
and A.L. [|^, The subject of the present paper is a generalization of the method developed in 
|2| to the general linear groups GL(N). The result will be a method to calculate the exponential 
of a quadratic matrix H, where only rational functions of the eigenvalues of H and the first N-1 
powers of H are involved. The key points are the Cayley-Hamilton theorem and the introduction 
of a multiplier m. 

The organization of this paper is as follows. First the method is shown in the low dimensional case 
SU(3) which is a group occurring quite often in physics. Then the general case of the expansion 
of the exponential map for a N x N matrix is presented. The method uses eignvalues 
The only real problem that remains is the determination of the eigenvalues of the matrix H. 
Throughout the paper we assume that the groups GL(N) are represented as N-dimensional ma- 
trices and that the eigenvalues of H are all different (cf. remark in appendix |B.3| ). 
Some possible applications of these results will be presented in a future paper. 



2 The exponential map for the group SU(3) 

The group SU(3) is used in several branches of physics. The best known application is the model 
of the strong interaction (see e.g. [^). For this reason and because its a good exercise to follow 
the steps of the general method we will show the exponential mapping of SU(3) in great detail. 
The calculations depend in some points on the fact that we consider a special group, i.e., the sum 
over the eigenvalues of the generator vanishes. There is no conceptional problem to extend the 
method to U(3). Like in the other cases (cf. fl], ^) a typical element U£ SU(3) can be written 
as exponential of the generator 

oo 

U = = h.H'' with H G su(3) . (1) 

?i=0 

The Cayley-Hamilton theorem and the iterated form in this case read 

ij3 ^ boH + Co and H^+' = aiH'^ + hH + a 

where the coefficients 6o and cq are functions of the eigenvalues of the eigenvalues x, y, z of H. 
They satisfy the recurrence relations 

Oj+i = bi , bi+i = ttibo + Ci , Ci+i = aj Co . (2) 
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Hence the coefficients satisfy 

fli+i = bo + aj_2 Co 

witli tlie first few values 

ao = 0, ai = bo, 02 = cq, as = 60^, «4 = ^0 co • 
Tlie explicit form of bo and cq can easily be derived from the secular equation 
= {X — x) {X — y) {X — z) = — {x + y + z) + {x y + X z + y z) X — X y z 



(3) 



ao = 



-60 



CO 



The leading coefficient oq vanishes since the generator H is traceless, i.e., x + y + z = 0. The 

+ t/^ + z'^). There are also some nice 
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2 X + y 



second coefficient can also be written as bo 
relations 

box + Co = x^ , boy + Co = y^ , boz + co = z^ . 
The idea now is to use the Cayley-Hamilton theorem for writing the sum (||) like 

U = h + H + \H'^ + Y.'^=o -^"^^ 

= h + H + + Y.n=0 ^ + hnH + C„) 

This form contains only sum over rational functions, there are no higher powers of the generator 
present anymore. The next step is now to find an analytic expression for the sums over the 
coefficients. 

A convenient form of the functions a„ , and c„ can be obtained if we introduce the multiplier 

m = {x — y) {x — z) {y — z) = {xP' {y — z) + y^ (z — x) + z^ {x — y)j . 
Then we get for the group element 



mU = m(h + H + + 



E 



m an 



mbn 
n=o(- + 3)!j 



H + 



mcn 



.S(- + 3)' 



It can easily be shown that the following form for the coefficients satisfy the recurrence relations 
m and (I) 



m an 


= {y - z) 


+ {z 


- x) 7/"+=^ 


+ {x 


- y) 


mbn 


= {y - z) 


+ {z 


- x) 2/"+^ 


+ (x 


- y) 


mCn 


= yz{y - 


z) + 


xz{z — x] 


yn+3 


+ xy{x - y) z^^^ 



The three sums are now 



E 

.n=0 



m an 
(n + 3) 



00 T 

^ mbn 
.nt^o(" + 3)! 



E 

.n=0 



mCn 

(n + 3)! 



{y - z) + {z - x) + {x - y) 
X {y — z) + y {z — x) + z {x — y) — 
y z {y — z) + x z {z — x) + xy {x — y) e 



m 
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Finally, we get the expansion of a SU(3) group element 



mU = [y z {y - 


- z) + X z {z — x) 


+ 


xy [x 


- y) el h 


+ {y 


- z) + y {z - x) 


+ 


z {x — 


y) e'] H 


+ Ky - 


z) + (z - x) + 


X 


- y) e'] 





3 The exponential map of GL(N) 

As we have seen in the cases of the groups SU(3) and SU(2,2) Q the exponential map can be 
written as sum over the first (N-1) powers of the generator H £ su(N). Where the coefficients 
are functions of the eigenvalues of H. In this section we generalize the results we have gotten on 
the low dimensional examples. It seems that there is a relatively easy concept of generalization. 
The desired result is an expansion of the exponential map of the form 

n=0 k=0 

where the coefficients are rational functions of the eigenvalues {A^; i = 1,2, ... , N} of the 
generator H. 



3.1 The secular equation 

The first step will be to take a look at the eigenvalues, some auxiliary functions, and their 
interrelations. 

Let us consider the secular equation of the matrix H 



N / N \ 

- A.) = - (T.CkX''-' 

1=1 \fc=0 / 



(6) 



where the coefficients are functions of the eigenvalues of H. In section 3.4 some coefficients 
are listed in their explicit form. For later convenience we introduce also the "truncated" version 
C(j);j of the coefficients C^, defined by 

N-l 

- A,) =: - E^W^A^-i-^ (7) 

j^i k=0 

Essentially the C(j)fc contain all terms of without A^. The connection between these coefficients 
can be seen easily via 
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N N 

n (A - \) = (A - Ai) n (A - A,) 

i=l 1=2 



(A - AO [\^-^ 
\^ + Ai C(i)7v-i 

N-l 

A^ — ^ Ck+i A 

A;=0 



N-l \ 

k=l / 
N-2 



fc=0 



,N-l-k 



Since the calculations above can be generalized to all eigenvalues we get the relations 
Ck = Cii)k - AjC(j)fc_i for A; = 1,2, ...,A^- 1 



Cat 



Ai C(j)Ar_i 



Let us define the multiplier m; i.e., the discriminant of the secular equation 

m := W (Ai - Aj) 

and the functions (see also (|30|)) 

m (Ai , . . . , Ai„i , Aj+i , . . . , \n) 



rrii 



n (^.- - ■ 



j<k 



(8) 



(9) 



(10) 



In what follows we mainly use the following form of m which can be obtained by expanding the 
Slater determinant (see (26) and ([29|)) 



N 



m 



i=l 



This formula can be generalized to 



N 



m 5ki 



J2 (-1)' m qi)fc-iAi^-^ for A;, / = 1, 2, . . . , iV . 



(11) 



(12) 



i=l 



For the proof see section 3.4. 



3.2 Recurrence relations 

The described method relies on the Cayley-Hamilton theorem which gives us the ability to write 
all powers H^^'^ for n G N in terms of the first N-l powers of H. The Cayley-Hamilton theorem 
for H G gl(N) reads 

N 

H"" = Y.'^kH^-'' (13) 
fe=i 
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The coefficients are the same as those in the secular equation @ and satisfy the recurrence 
relations ( p^ derived below. For the special groups, i.e., det g = 1 for g G SL(N) the first 
coefficient vanishes since the sum over the eigenvalues is zero. 



Multiplication of (13) with and using (|l^) again gives the iterated form 

N 
k=l 

Multiplying once more with H gives 

+ ... + (C+i + crc„)i?^-" + ... + 

+ (Cat + Ci Cn-i) H + C^Cn 

N 

_L ^ ^n+l jjN-k 
k=l 



(14) 



and hence we get the recurrence relations 



N-l — '^N + '-^l '-'N-l 



Cn+1 /^n I /^n 
2 — O3 + U]^ U2 

Ul Gat . 



(15) 



If we successively plug in the in the recurrence relation of Cf we get a formula which contains 
only terms with Cf and the coefficients of the original Cayley-Hamilton equation (O) 



CI 



N 

J2 Cr^^~^' Cj for n > iV - 1 

n 

^ C""-'' Ci+j + m C7„+2 for n < iV - 1 
I 3=0 



(16) 



For the other coefficients C^+^ {k = 1,2,...,N) we get analogous formulae 



N-k 



Cl'-^ Ck+j for n > N-k 



j=0 



(17) 



^ C" Cfc+j + m Ck+n+i for n < N - k 
I i=o 

The coefficients of the secular equation have the explicit form 

(11) " 



m Ck 



j=l 



TV 



1=1 

N 



N-l 



C, 



(i)fc 



2=1 



i=l 



: for fc 7^ 
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where we applied Eq.(34) to the second term in the last equation. We will need this form as first 
values in the proof of Eq.(21) 



N 



m Ck 



^ (-1)^ nii Xf for A; = 1, 2, . . . , iV 



i=l 



(18) 



From the SU(3) and SU(2,2) cases one may assume that the recurrence relation ( [iq ) has the 
solution 



N 



N+n 



i=l 



Proof: 

The proof of Eq.(p!9|) is done by induction over n. 
The first coefficient (n = 0) is given by 



N 



mCi = ^(-iy+^m,A, 



N 



i=l 

which is easy to prove if one writes Ci in the form 

N 



Ci 



Xj = Xi + = K + C(i)i . 



For the product m Ci we take m in the form of Eq.(29) we get 



N 



m 



1=1 

N N 



N-l 



(19) 



(20) 



The last equation holds since the exponent of Aj should be — 2 in order to yield a non- vanishing 
sum (see Eq. (|3^) ) . 



First we treat the case of n > A^. The next step is to assume the validity of ( |19D for n and to 
show that then it follows also for n + 1 



^n+l (16) \ - „n+l-j ^ © 



N N 

j=i i=i 



X N N 

j=i i=i 



6 



N 



N 



AT 

E 

i=i 



y+l 

i=l \j=l 
N ( N+1 N 



i=l 



\i=2 



i=l 
N 



\ n+1 t~i \ Af+n+l 



i+1 ^ . \ ,JV+ra+l 



In the case of n < A'^ — 1 there is an additional term 

n 

AT 



Af AT 
i=l i=\ 



Y (-1)'+' Af +"+1 + Y (-1)'^' Af + m C, 



n+2 



mC„4 



(18) 



AT 



m.A 



i+l ^ . -v A^+n+l 



i=l 



The coefficients C|? can be written in the form 



mCl 



N 

Y i-^y m q.)fc-iA.^+" 



for k = 1,...,N 



□ 



(21) 



Proof 



The proof is analogous to the one of m Cf but uses the explicit form (^) of these coefficients. 



j=0 

N-k N 



j=0 i=l 



Y E (C(i)fc+j - Aj C(i)A;+j_l] A, 

j=0 i=l 



N+n-j 



N-k N 

j=0 i=l 

N~k N 



i=l 



If we now shift the summation index j in the second sum most of the terms cancel with those 
of the first sum. The remaining term j = N — k in the first sum contains C(j)jv = 0) hence, 
vanishes also. Therefore, only the third sum remains what proofs the assumed form (^Tj) of C^. 
Again there are the cases n < N — k which need to be treated separately 

N-k 

™ = E C-i"-^' Ck+j + m Ck+n+i = ...= 

j=0 

N N 

= E (-1)' ^^ Ci^k-l Af + E (-1)'^' C^r)k+n Af + ^ 
i=l i=l 



= -mCfc+„+i (18) 

N 

E (-1)^ m C^i)k-i Af . 

i=l 



□ 



3.3 The exponential map 

The expansion of a group element g G G with generator H € gl(N) can now be written like 



E 

n=0 



n=0 



1 



n=0 



n 



{N + n)l 
E C'k H^-'' 



n=0 
/ A' 



^k=l 



Using the multiplier m we get 

E ^ + E 



E 



1 



n=0 k=l Ln=0 ^ ' 

We can now treat the sums for different k separately 



N-k 



E 

n=0 



{N + ny. 



mC2 



E 



1 



N 



u^+^yti 



E (-1)' m C7(,)fc_i A, 



N+n 



(22) 



(23) 
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N 



1=1 n=0 ^ ' 

N / N-1 ^ „\ 

E (-1)^ - E ^ 

i=l \ n=0 / 

TV N-1 , TV 

E (-1)' + E ;;t E (-1)'^' 
1=1 



Af-l N 

$Z „ I E/ 

n=0 i=l 



The last equation relies on Eq. (^3|) and (|^) . The terms — ^jy'^^^, cancel the first sum in Eq. (^) . 
The final result turns out to be 



m e 



^ = (-1)^ detF (-1)* rrii ^] In + (e %^-2 e"'") ^ 



v,i=l 



\i=l 



+ . . . + E (-1)' ^w'^ e^') + 



^+1 m,- i/^-i 



or in closed form 




(24) 



which reads in terms of the adjoints of the Slater determinant 



N-1 / N 



m e 



H 



E E ^(On 



n=0 \i=l 



(25) 



3.4 The Slater determinant 



One crucial ingredient of the method is the usage of a multiplier m, defined in Eq.(P). From low 
dimensional examples one may assume the forms ( pTj ) and (|l^) of m. The general proofs can be 
done by writing m as Slater determinant. The Slater determinant is defined as (cf. P) 



1 


1 


1 


\n 


Aat-i . 


• Ai 


^N 




• A? 


N-1 
'N 


\N-1 
^N-1 ■ 


• Af- 



n 



A,) 



(26) 
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We can now use the Laplacian method of expanding the Slater determinant 

n n 
det A = ^ ttj-j An = ^ ^ dij A-ij 



O-ij ^ij 



(27) 



i=l 



where the so-called adjoints Aij are the subdeterminants of aij multiplied by the sign factor 

{-iy+^. 

It is also well known that the Laplace expansion with "wrong" adjoints gives zero 



= ^ Qij Aij for I ^ i . 



(28) 



The Laplacian method applied with respect to the last row gives then the expansion (11) 



N 



N 



m 



N-l 



A 



{N+l-i)N 



1=1 

N 
i=l 



i=l 



N-l 



where m,- are the subdeterminants of m 



rrii 



n (^fc - ^i) • 



(29) 



(30) 



We get 



m 



(7) 



mi Jl (Afc - \i) W {\i - Xj) 

k<i i<j 

i-iy-' rrii n (Ai - Aj) 

AT-l 

(-l)-i ^ (- C7(,)„ A,^-^- 



for i = l,2,...,N 



n=0 



Af-1 



N-l 



E (-1)' ^1 



i ^{i)N-l-n A" 



E "^(«)« Aj' 



(31) 



n=0 



n=0 



what proves Eq.(^). Also Eq.(34) is proven since if the exponent of Aj is not n the sum vanishes 
because it is an expansion with the "wrong" adjoints ^(j)^- Hence we get an explicit expression 
for the adjoints 



A 



{i)n 



{-ly mi C(j)7V-n-l • 



(32) 



We can also expand m with respect to the (n+l)-th line and then use Eq.(32) 



N 



m 



J2 Ar for n = 0,l,...,iV-l 



i=l 
N 



E (-1)' 



i '^(i)N-l-n Aj" 



(33) 



i=l 
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Writing the Laplacian expansion with "wrong" adjoints leads to 

N 

= for k,n = 0,1,... ,N -I k^n 



i=l 

N 



^ (-1)* rUi C(j)Ar_i_fc A" 



i=l 



Ergo (cf. Eq.(|l2D) 



N 



m 5ki 



J2 (-1)' m for k,l = 1,2,...,N 



(34) 



(35) 



Appendix 

A Some details 

This section contains exphcit forms of some coefficients and some proofs. Almost all of the 
equations hold in the general case, but those which hold only in the case of the special groups, 

g 

i.e., vanishing sum of eigenvalues, are denoted by the sign = 
For the coefficients of the secular equation we get, e.g., 

N 

Cj, = (-1)^+1 = (-1)'^"'' deti/ 

1=1 

Cn^i = (-l)^En^. = (-l)^E^ (36) 

i=l jyti 1=1 * 

1 ^ 

C2 = (-1) E ^ 2 ^ ^'^^ 

i<j k=l 

N 

C^i = E = ' Co = - I. 

i=l 

Some "truncated" coefficients 

%iv-i = (-1)^"^ n A, = (-1)^ ^ 
C{i)N-2 = (-1)^"^ E n 



C(i)2 = (-1) E ' '^(01 = E 

C(i)0 = - 1 , C'(j)Ar = . 
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B Additional checks 



B.l One-dimensional subgroups 

One-dimensional subgroups (cf. |^]) of GL(N) can be generated by 

{e*^ ; H e gl(N) , t G r} . (37) 

From the known expansion of we can derive the expansion of e* ^ by multiplying the occurring 
expressions with an appropriate factor. Obviously, the eigenvalues of the t-dependent generator 
t H are t Aj if the Aj are the eigenvalues of H. Therefore, we need to make the replacements 

Aj — > t Aj 

Cfc — > Ck C{i)k — ^ t'' 

The expansion ( |25| ) reads now 

N-l / N \ 



n=0 \j=l / 
N-l / N \ 



n=0 \i=l / 

or 

N-l / N \ 

^^'"^ = E E(-l)^%n-imie*^ i/^-". (38) 

n=0 \j=l / 

Differentiation of the r.h.s of Eq.(|38|) and setting t = gives the derivation of the unit element 

N-l N 

E ( E (-1)' C^i)n-i m A. ) H""-^ = mH. 

n=0 i=l 



Since this result coincides with the one we get by differentiating the l.h.s it is an additional proof 
of the expansion (24). 



B.2 Eigenvalues 

It is easy to demonstrate that Eq.(p5|) gives also the right connection between the eigenvalues 
of the generator H and the ones of the corresponding group element g = . Let be the 
eigenvectors of H with eigenvalues Aj 

H Xi = \i Xi for z = 1, 2, . . . , . 

For the powers of H we get 

i?" Xi = A " Xi for n € N . 
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Plugging this in Eq.(p4|) yields 



N-l / N \ N-l 

me^ Xj = ( E ^{i)n ) H'^Xj = ^ ( A,") Xj 

n=0 

:i = m e^^ Xi . 



n=0 


\i=l 


N 


N-l 


E( 


E 


i=l 


n=0 




m 



Therefore, we get the desired result 

which again confirms the expansion (p^). 
B.3 Remark 

In the cases where some of the eigenvalues coincide, the multipher m will be zero. But in 
these cases m can be chosen in a simpler fashion so that there occur only non-vanishing factors. 
Essentially all factors which will become zero can be canceled out in Eg . (|2^) . 

Acknowledgement 

The author would like to thank Dr. K-P. Marzlin for the discussions of some points. 

References 

[1] A.O. Barut, J.R. Zeni, and A. Laufer. The exponential map for the conformal group 0(2,4)- 
J.Phys.A:Math.Gen., 27:5239-5250, 1994. 

[2] A.O. Barut, J.R. Zeni, and A. Laufer. The exponential map for the unitary group SU(2,2). 
J.Phys.A:Math.Gen., 27:6799-6805, 1994. 

[3] Ta-Pei Cheng and Ling-Fong Li. Gauge theory of elementary particle physics. Oxford Uni- 
versity Press, 1984. 

[4] Sergio Fubini. Vertex operators and Quantum Hall Effect. Mod. Phys. Lett., A6: 347, 1991. 

[5] Sigurdur Helgason. Differential Geometry, Lie Groups, and Symmetric Spaces. Academic 
Press, 1978. 

[6] Claude Itzykson and Jean-Michel Drouffe. Statistical Field Theory II. Cambridge University 
Press, 1989. 

[7] J.Ricardo Zeni and Waldyr A. Rodrigues. Hadronic Journal, 13:317, 1990. 

[8] J.Ricardo Zeni and Waldyr A. Rodrigues. A thoughtful study of Lorentz transformations by 
Clifford algebras. International Journal of Modern Physics A, 7(8): 1793-1817, 1992. 



13 



